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The effect of medium-induced parton energy loss on jet fragmentation is studied in high-energy 
heavy-ion collisions. It is shown that an effective jet fragmentation function can be extracted from 
the inclusive pt spectrum of charged particles in the opposite direction of a tagged direct photon 
with a fixed transverse energy. We study the modification of the effective jet fragmentation function 
due to parton energy loss in AA as compared to pp collisions, including Et smearing from initial 
state radiations for the photon-tagged jets. The effective fragmentation function at z = pr/Ef ~ 1 
in pA collisions is shown to be sensitive to the additional Et smearing due to initial multiple parton 
scatterings whose effect must be subtracted out in AA collisions in order to extract the effective 
parton energy loss. Jet quenching in deeply inelastic lepton-nucleus scatterings as a measure of the 
parton energy loss in cold nuclear matter is also discussed. We also comment on the experimental 
feasibilities of the proposed study at the RHIC and LHC energies and some alternative measurements 
such as using Z° as a tag at the LHC energy. 


I. INTRODUCTION 

Hard processes are considered good tools to study ul- 
trarelativistic heavy-ion collisions because they happen 
early in the reaction processes and thus can probe the 
early stage of the evolution of a dense system, during 
which a quark-gluon plasma (QGP) could exist for a 
short period of time. Among the proposed hard probes, 
large transverse momentum jets or partons are especially 
useful because they interact strongly with the medium. 
For example, an enhanced acoplanarity and energy im¬ 
balance of two back-to-back jets JtJ due to multiple scat¬ 
terings, jet quenching due to the medium-induced ra¬ 
diative energy loss of a high-energy parton propagating 
through a dense medium [Q can provide important infor¬ 
mation on the properties of the medium and interaction 
processes that may lead to partial thermalization of the 
produced parton system. The medium-induced radia¬ 
tive energy loss of a fast parton traversing a dense QCD 
medium is also interesting by itself because it illustrates 
the importance of quantum interference effects in QCD. 
As recent studies have demonstrated [HU, it is very im¬ 
portant to take into account the destructive interference 
among many different radiation amplitudes induced by 
multiple scatterings in the calculation of the final ra¬ 
diation spectrum. The so-called Landau-Pomeranchuk- 
Midgal effect || can lead to very interesting, and some¬ 
times nonintuitive results for the radiative energy loss 
of a fast parton inside a QCD medium. Recently, Baier, 
Dokshitzer, Mueller, Peigne and Schiff (BDMPS) showed 
|ij that the energy loss per unit distance, dE/dx , grows 


linearly with the total length of the medium, L, which 
in turn can be related to the total transverse momentum 
broadening squared, Afc^, of the parton from multiple 
scatterings. It turns out that both quantities, dE/dx and 
are related to the parton density of the medium 
that the parton is traveling through. One can there¬ 
fore determine the parton density of the produced dense 
matter by measuring the energy loss of a fast parton in 
high-energy heavy-ion collisions. 

Unlike in the QED case, where one can measure di¬ 
rectly the radiative photon spectrum and thus the energy 
loss of a fast electron, one cannot measure directly the 
energy loss of a fast parton in QCD. Since a parton is ex¬ 
perimentally associated with a jet, a cluster of hadrons 
in a finite region of the phase space, an identified jet 
can contain particles both from the fragmentation of the 
leading parton and from the radiated partons. If we ne¬ 
glect the kp broadening effect, the total hadronic energy 
contained in a jet should not change even if the leading 
parton suffers radiative energy loss. However, significant 
changes could happen to particle distributions inside the 
jet, or the fragmentation function and jet profile, due to 
the induced radiation of the leading parton. Therefore, 
one can measure the radiative energy loss indirectly via 
the modification of the jet fragmentation function and 
jet profile. 

A jet fragmentation function is defined as the particle 
distribution in the fractional energy. In order to measure 
the fragmentation function one has to first determine the 
initial energy of the fragmenting parton either through 
other measured kinematic variables as in e + e _ and e~p 


1 


or calorimetric measurements as in pp and pp collisions. 
However, because of the large value of dET/dyd<j> and 
its fluctuation in high-energy heavy-ion collisions || , the 
conventional calorimetric study of a jet cannot determine 
the jet energy to such an accuracy as required to deter¬ 
mine the parton energy loss. In search for an alterna¬ 
tive measurement, Wang and Gyulassy JtJ proposed that 
single-particle pt spectrum can be used to study the ef¬ 
fect of parton energy loss, since the suppression of large 
Et partons naturally leads to the suppression of large 
Pt particles. Because the particle spectrum at large py 
is the convolution of jet cross sections and fragmentation 
functions, particles with a fixed pt can come from the 
fragmentation of partons of different initial energies with 
some average value ( E 1 ^ t ). In this case, the suppression 
of the pt spectrum due to parton energy loss is then re¬ 
lated to the modification of jet fragmentation function 
at an averaged (z) = Pt/(E 3 ^ t ). Since (Ey^) is approx¬ 
imately proportional to pt, by varying pt one can then 
study the energy dependence of the modification of jet 
fragmentation functions at a fixed (z) ||. 

In order to study the modification of the whole jet frag¬ 
mentation function due to parton energy loss in the full z 
range, we and Sarcevic || proposed to measure the par¬ 
ticle pt distribution in the opposite transverse direction 
of a tagged direct photon. Since a direct photon in the 
central rapidity region (y = 0 ) is always accompanied by 
a jet in the opposite transverse direction with roughly 
equal transverse energy, the pt distribution of particles 
in that direction is directly related to the jet fragmen¬ 
tation function with known initial energy, Ey 4 ~ Ey. 
In such 7 + jet events, the background due to particle 
production from the rest of the system was estimated to 
be well below the pt spectrum from jet fragmentation at 
moderate large py. Therefore, one can easily extract the 
fragmentation function from the experimental data with¬ 
out much statistical errors introduced by the subtraction 
of the background. By comparing the extracted jet frag¬ 
mentation function in AA to that in pp collisions, one 
can then measure the modification of the fragmentation 
function and determine the parton energy loss. 

Because of the complexity of the problem, it will be 
helpful for us to first discuss all possible relevant pro¬ 
cesses which might have some effects on the study of 
parton energy loss in 7 + jet events in high-energy AA 
collisions. One can order the processes in a chronological 
order with respect to the hard process of direct photon 
production. 

( 1 ) As the two nuclei approach and pass through each 
other, the two participating beam partons which later 
produce the direct photon will suffer initial state inter¬ 
actions with other oncoming nucleons. The participating 
beam partons will then suffer radiative energy loss and 
acquire transverse momentum kicks because of these soft 
interactions. The initial state interactions will also cause 
the shadowing of the parton distributions inside the nu¬ 


clei. These initial state effects will certainly affect the 
production rate of direct photons at a given transverse 
energy Ej. However, they will not influence the propa¬ 
gation and fragmentation of the accompanying produced 
jet parton in events triggered with a direct photon with 
a fixed E y. 

(2) After the hard process in which a direct photon 
and a jet parton are produced, the jet parton can also 
scatter from the beam nucleons within a tube of a trans¬ 
verse size at most 1 fm in the rest part of the colliding 
nuclei which has not passed through. Since the colliding 
nuclei pass through each within a very short period of 
time, t ~ 1 fm/c (this is the spatial size of wee-partons, 
while the valence partons have a Lorentz contracted size 
of RATn 2 N / 2 s), the produced jet parton in central rapid¬ 
ity region will not have time to interact with other beam 
nucleons outside the tube. Since we are only interested 
in jet partons in the central rapidity region, these scat¬ 
terings will not cause the jet partons to lose transverse 
energy. Rather, together with the initial state interac¬ 
tions, they will change the final transverse momentum 
of the jet parton, resulting in an Et broadening in addi¬ 
tion to the Ey-smearing caused by initial state radiations 
associated with the hard process. The effects of this Et 
broadening should also exist and can be studied indepen¬ 
dently in pA collisions. 

(3) In the triggered events, there are many other pro¬ 
cesses which can also produce hard or semihard partons 
0 and thus form a dense medium in the central rapid¬ 
ity region. The photon-tagged jet parton will then in¬ 
teract with these partons during its propagation through 
the dense medium. We call these interactions as final 
state interactions. The induced radiative energy loss 
and transverse momentum broadening, referred to as fey- 
broadening in this paper, are the focus of our study. 

In this paper, we will study in detail the effect of par- 
ton energy loss on the jet fragmentation function as ex¬ 
tracted from the pt spectrum in the opposite direction of 
a triggered direct photon. In particular, we will take into 
account the Et smearing of the jet due to initial state 
radiations associated with the 7 + jet processes. We will 
show that the particle spectrum from the jet fragmen¬ 
tation at pt ~ Ey is very sensitive to the Et broad¬ 
ening from initial and final state scatterings with beam 
partons. One can then use our proposed measurement 
to determine the Et broadening in pA collisions. This 
small but finite effect must then be subtracted out when 
one determines the medium-induced parton energy loss 
in A A collisions. We will also investigate the sensitivity 
of the modification of the fragmentation function to the 
energy and A dependence of the parton energy loss. The 
change of the profile function in the azimuthal angle due 
to the fey-broadening of the parton from multiple scat¬ 
terings inside the medium will also be discussed. Here fey 
is the parton transverse momentum with respect to the 
original jet direction, which can be related to the parton 
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energy loss according to BDMPS study [Q. Finally, we 
will discuss the experimental feasibility of the proposed 
study and alternative measurements using Z° particles 
as a tag. We will also discuss how similar measurements 
can be made in deeply inelastic lepton-nucleus scatter¬ 
ings, from which one can determine the energy loss of a 
fast parton passing through a cold nuclear matter. 

II. MODIFIED FRAGMENTATION FUNCTIONS 

The fragmentation functions of partons hadronizing in 
the vacuum have been studied extensively in e + e _ , ep 
and pp collisions (jllj. These functions describe parti¬ 
cle distributions in the fractional energy, z = Eh/Ej et , 
in the direction of a jet. Similar to parton distributions 
inside hadrons, the fragmentation functions are also non- 
perturbative in nature. However, parton cascades during 
the early stage of the fragmentation can be described 
by perturbative QCD. The measured dependence of the 
fragmentation functions on the momentum scale is shown 
to satisfy the QCD evolution equations very well. We 
will use the parametrizations of the most recent analysis 
|n| in both z and Q 2 dependence for jet fragmentation 
functions Z?°^ a (z, Q 2 ) to describe the fragmentation of a 
parton (a) into hadrons ( h ) in the vacuum. 

The fragmentation of a parton inside a medium is dif¬ 
ferent from that in the vacuum, because of its final state 
interactions with the medium and the associated radia¬ 
tions. Such interactions and medium-induced radiations 
will cause the deflection and energy loss of the propagat¬ 
ing parton which in effect will modify the fragmentation 
functions from their corresponding forms in the vacuum. 
In principle, one could study the modification of jet frag¬ 
mentation functions in perturbative QCD in which in¬ 
duced radiation of a propagating parton in a medium 
and Landau-Pomeranchuk-Migdal interference effect can 
be dynamically taken into account. However, for the pur¬ 
pose of our current study, we can use a phenomenological 
model to describe the modification of the jet fragmenta¬ 
tion function due to an effective energy loss dE/dx of the 
parton. Such an approach is useful and possibly neces¬ 
sary for experimental studies of the parton energy loss 
and multiple final state scatterings. 

In this phenomenological model we assume that the 
size and life time of the system is small compared to the 
hadronization time of a fast parton. In the case of a 
QGP, a parton cannot hadronize inside the deconfined 
phase. A fast parton will hadronize outside the system 
and the fragmentation can be described as in e + e _ colli¬ 
sions, however, with reduced parton energy. The interac¬ 
tion of a parton a with the medium can be characterized 
by the mean-free-path A a of parton scatterings, the ra¬ 
diative energy loss per scattering e a and the transverse 
momentum broadening squared AAfy. The energy loss 


per unit distance is thus dE a /dx = e a /X a which in prin¬ 
ciple depends on Afcf. and A a . We assume that the prob¬ 
ability for a parton to scatter n times within distance L 
is given by a Poisson distribution, 

a(n,£) = fGMb-w.. (i) 

n! 

We also assume that the mean-free-path of a gluon is half 
that of a quark, and the energy loss dE/dx is twice that 
of a quark. The emitted gluons, each carring energy e a 
on the average, are assumed to hadronize also accord¬ 
ing to the fragmentation function. For simplification, we 
will neglect the energy fluctuation given by the radiation 
spectrum for the emitted gluons. We assume the momen¬ 
tum scale in the fragmentation function for the emitted 
gluons to be set by the minimum scale Q q = 2.0 GeV 2 . 
Since the emitted gluons only produce hadrons with very 
small fractional energy, the final modified fragmentation 
functions in the moderately large z region are not very 
sensitive to the actual radiation spectrum and the mo¬ 
mentum scale dependence of the fragmentation functions 
for the emitted gluons. In this paper, we will also neglect 
possible final state interactions between hadrons from 
parton fragmentation and the hadronic environment at 
the late stage of the evolution of the whole system. How¬ 
ever, it is important for future investigations to estimate 
the influence of pure hadron scatterings on the final ob¬ 
served jet fragmentation functions. 

We will consider parton fragmentation in the central 
rapidity region of high-energy heavy-ion collisions. In 
this case, we only need to study partons with initial trans¬ 
verse energy Et and traveling in the transverse direction 
in a cylindrical system. With the above assumptions, the 
modified fragmentation functions for a parton traveling 
a total distance L can be approximated as Q, 

i N a 

D h/ a(z,L,Q 2 ) = — J2Pa(n, L)—D^ /a (z“, Q 2 ) 
n=0 

+ (n a )^D° h/g (z' a ,Q 2 ), (2) 

where z“ = z/(l - ne a /E T ), z' a = zE T /e a and C% = 
Y^,n=o Pa{n)- D°/ q (z, Q 2 ) are the jet fragmentation func¬ 
tions in the vacuum which we take the parametrized 
form in Ref. ]l^]. We limit the number of inelastic scat¬ 
terings to N = Et /e a to conserve momentum. One 
can check that the above modified fragmentation func¬ 
tions satisfy the momentum sum rule by construction, 
E h f zD h/a{z,L,Q 2 )dz = J2 h f zD h/a( z ’Q 2 ) dz = !• 
For large values of N, the average number of scatter¬ 
ings within a distance L is approximately (n a ) ~ L/X a . 
The first term in the above equation corresponds to the 
fragmentation of the leading partons with reduced en¬ 
ergy Et — ne a after n inelastic scatterings. This term 
normally dominates for leading particles in the moderate 
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and large z region. Since the fragmentation functions 
D^/ a { z ,Q 2 ) generally decrease with z, especially quite 
rapidly at moderate and large z region, the reduction in 
energy will lead to the suppression of leading particles or 
the decrease of the fragmentation functions in this region 
as compared to the case in vacuum. The second term in 
the above equation comes from the emitted gluons each 
having energy e a on the average. This term is generally 
significant only in the small z region and it increases the 
effective fragmentation functions in the small z region, or 
enhances soft particle production. We should note that 
our assumptions on the hadronization of the emitted glu¬ 
ons are too schematic to give a quantitative description 
of the physics involved in that small z region. 

For a given parton energy Et and the total distance 
L, the above effective fragmentation functions depend 
on only two parameters, the mean-free-path A a and en¬ 
ergy loss per scattering e a . As demonstrated in Ref. 
contributions from the leading partons who have suf¬ 
fered at least one inelastic scattering is completely sup¬ 
pressed for 2 values close to 1. The remaining contri¬ 
bution comes from those partons which escape the sys¬ 
tem without a single inelastic scattering, with a prob¬ 
ability exp(— L/X a ), which depends on A a but is inde¬ 
pendent of the parton energy Et and the parton energy 
loss dE a /dx. On the other hand, in the intermediate z 
region, particles from the fragmentation of the leading 
partons with reduced energy dominates. The suppres¬ 
sion of the fragmentation functions is controlled by the 
total energy loss, (A Et) = ( n a )e a = LdE a /dx, which de¬ 
pends only on dE a /dx. One, therefore, could determine 
in principle these two parameters, A a and dE a /dx , simul¬ 
taneously from the measured suppression of the effective 
fragmentation functions, for fixed Et and L. However, 
as we will see in the next section, the complication of 
not knowing the jet energy precisely will render such ar¬ 
guments unrealistic. In certain cases, one has to resort 
to a model-dependent global fitting of the modification 
of the fragmentation functions in order to determine the 
mean-free-path and parton energy loss. 


III. THE INCLUSIVE FRAGMENTATION 
FUNCTION OF PHOTON-TAGGED JETS 

As we have emphasized in the Introduction, the most 
important point in the study of parton energy loss 
through the measurement of the modification of the par- 
ton fragmentation functions is the determination of the 
initial parton or jet energy. However, the direct mea¬ 
surement of a jet energy to the accuracy as required to 
determine an energy loss of a few GeV is unfeasible due to 
the large background and its fluctuation in high-energy 
heavy-ion collisions. To overcome this difficulty, it was 
proposed Q that direct photons in heavy-ion collisions 
can be used to tag the energy of jets which always accom¬ 


pany the direct photons. Because of initial state radia¬ 
tions associated with the production of a direct photon, 
the accompanying jet is not always exactly in the oppo¬ 
site direction of the photon and its transverse energy also 
differs from collision to collision, though the averaged jet 
energy is well approximated by the energy of the trig¬ 
gered photon. In Ref. the variation of jet energy was 
not considered in the study of the effective inclusive jet 
fragmentation function in 7 + jet events and the modifi¬ 
cation due to parton energy loss. In this paper, however, 
we would like to explore the effect of Et smearing due 
to initial state radiations. We will see that such a smear¬ 
ing complicates the simple procedure to determine the 
parton energy loss and the mean-free-path of final state 
scatterings as outlined in the previous study ||. 

Let us consider events which have a direct photon with 
fixed transverse energy E/ in the central rapidity region, 
|y 7 | < Ay/2, Ay = 1. Given the jet fragmentation func¬ 
tions D/,(z), with z the fractions of momenta of the jet 
carried by hadrons, one can calculate the differential pt 
distribution of hadrons, averaged over the kinematical 
region (Ay, A(j>), from the fragmentation of a photon- 
tagged jet in pp collisions, 


dN& 


-h d 


d&pp 


-h J 


da/ p 


dyd 2 pT 
± 


-h* 


dy 1 dE/.dyd 2 pT 
D° h/a (PTlE “) 

x — - -—- 

PtEZh 


da/p / dy-y dE/ dy 7 dE/ dyd 2 pr 
= / dE/dy a d</>, 


da/ p a 


1 dE/dy^dE/dyad/a 


l lT-ir7fo(ya-y,<t>a-(l>) (3) 

I (Ay,Atp) Ay A(j) 

where the summation is over jet (a) and hadron ( h ) 
species and fo(y, />), assumed to be the same for all 
hadron species, is the normalized hadron intrinsic pro¬ 
file around the parton axis. If the azimuthal angle of 
the photon is 0 7 and </> 7 = </> 7 + n, the restricted kine¬ 
matical region for the selected hadrons is defined as 
(Ay, A<j>) = (|y| < Ay/2,\<j> - 0 7 | < A<j>/2). One could 
also use more complicated geometry, such as a circle with 
a given radius, for the phase space restriction to define 
a jet. The inclusive differential cross section for direct 
photon production is 


do/p 

dy 1 dE/ 


E 


dE/dy a d<f) a 


da/p a 


dE/ dy-y dE/ dy a d(j) a 


( 4 ) 


We now define the F/r-smearing function, g pp (E/.E/), 
and parton correlation function, /jet(y<j, </>a), as 


9pp(E T i E/)fj et (y a , (j> a ) — 


da^~ a 

UUpp 


da vp dE/dy^dE/dyad/a 

dy-, dE J 


( 5 ) 

In a perturbative calculation to the lowest order in 
a s , which was used in our earlier study ||, the Fir- 
smearing function and parton correlation function are 
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simply two ^-functions, g pp (E^, Ej,) = S(E^ — Ej,), 
f jet (y a, (fra) oc 5{<fr a — 0 7 ), due to momentum conserva¬ 
tion. The intrinsic hadron profile function fo(y,(fr) in 
this leading order calculation should be the measured jet 
profile. In calculations beyond the leading order, the pho¬ 
ton and jet parton have a finite imbalance in transverse 
momentum due to the initial state radiations. The final 
state radiations also contribute to the measured jet pro¬ 
file f(y, (fr) which should be the convolution of the parton 
correlation function /j e t(z/a, (fra) from perturbative calcu¬ 
lations to a given order and the intrinsic hadron profile 

fo(y,<fr), 

f{y, (fr)= dy a d(frafjet(ya, (fra)fo(y a -y,(fra-(fr) • ( 6 ) 


Note that the differential cross section da 7 ~ a and the 
fragmentation function D^ a (z,Q 2 ) in Eq. (||) both de¬ 
pend on the factorization scheme and the associated 
scale. So are the .Bp-smearing function and parton cor¬ 
relation function in Eq. (||). Unlike the collinear di¬ 
vergences in the next-to-leading order jet cross sections, 
which are cancelled via the definition of a jet with a finite 
size, the collinear divergences in Eq. (||) are subtracted 
out via the definitions of parton distributions and frag¬ 
mentation functions. Therefore, the differential cross sec¬ 
tion der 7_a beyond the leading order in Eq. (|j) is not the 
same as the 7 — jet cross section in which a jet is de¬ 
fined via the transverse energy within a finite region in 
phase space 0. The two cross sections can be related 
via some divergence-free physical observables, e.g., total 
hadronic energy within the acceptance (Ay,A(fr), which 
can be computed from Eq. (|S|). 

With the above definitions of the Bp-smearing function 
and jet profile function, we can rewrite Eq. ((3j) as 


dN;~ h± 

dyd 2 pr 


^V 0 (Bj) J dE T g pp (E T ,E^) 

a,h 


D h/a(PT/E T ) C(Ay,A<fr) 
PtEt AyAcfr 


(7) 


where C(Ay, A<fr) = f M < Ay/2 dy ^_^,< A0/2 d(frf{ Vl (fr- 
</> 7 ) can be considered as an overall acceptance factor for 
finding the jet fragments in the given kinematic range. 

We will approximate the fractional production cross 
section, r 0 (Bj), of a-type jet associated with the direct 
photon, by the lowest order calculation, 


r a (El) 

do 2 

dy-fdE t 


daZ/dy^dE^ _ ^ 
do 1 /dy^dEfrp 
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d%bfb/p(.Xb) fc/p(.Xc )" 


be Jx bmin 


X b X c 


dcT 


2xb — xre v -y dt 


( 6 c-> 7 + 0 ), 


( 8 ) 


(9) 


where x c = XbXT£~ v ' 1 / {2xb — xpe^), Xb m in = xpE 1 ' 1 / {2 — 
xt&~ V i ), and Xt = 2Ep/y/s . The parton distribu¬ 
tions in a proton, f a / P (x ), will be given by the MSRD—/ 
parametrization [ fj~7| . In our following calculations for 
AA collisions, we will use the impact-parameter averaged 
parton distributions per nucleon in a nucleus (A, Z), 


fa/A(x) = S a/A (x) 


z z 

~^fa/p{x) + (1 J )fa/n{x) 


A' 


, (10) 


where S a / A (x) is the parton nuclear shadowing factor 
which we will take the HIJING parametrization jl3| . We 
have explicitly taken into account the isospin of the nu¬ 
cleus by considering the parton distributions of a neutron 
which are obtained from that of a proton by isospin sym¬ 
metry. 



FIG. 1. The normalized parton correlation from HIJING 
simulations in rapidity y and azimuthal angle <f> with respect 
to the opposite direction of a tagged photon with = 10 
GeV in pp collisions at y^s = 200 GeV 

To simulate higher order effects, event generators such 
as PYTHIA |l5[, which was used in HIJING program 
p[ , normally use parton shower model. In this model, 
one introduces a cut-off /io for the parton virtuality in the 
chain of parton shower, thus avoiding both infrared and 
collinear singularities. In addition, initial and final state 
radiations are treated separately and the interference be¬ 
tween them is also neglected. Shown in Fig. [t], is the par- 
ton correlation function in rapidity y and azimuthal an¬ 
gle efr with respect to the opposite direction of a triggered 
photon with E J. = 10 GeV in HIJING |li| simulations 
of pp collisions at ^/s = 200 GeV. In the simulations, 
the final state radiations are switched off so that we can 
study the effect of momentum imbalance due to initial 
state radiations. We can see that most of the jet partons 
fall into the kinematic region, (\y\ < 1, |</> — 0 7 | < 0.5). 
The jet profile, which is the convolution of the parton 
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correlation function and hadron intrinsic profile around 
the parton axis, has a similar shape with a slightly larger 
width according to both our simulations and experimen¬ 
tal measurements in high-energy pp collisions |h]]. For 
calculations throughout this paper, we will use A y = 1 
and At/) = 2. We find the acceptance factor defined in 
Eq. (|^), C(Ay, At/>) ss 0.5, independent of both the col¬ 
liding energy and the photon energy EJ, using HIJING 
[|l3| Monte Carlo simulations. 



E t (GeV) 

FIG. 2. The i?T-smearing function for the photon-tagged 
parton jets with E J = 10, 15 GeV, from HIJING simulations 
of pp collisions at s/s = 200 GeV. The averaged value of Et 
is (Et) = 8.08, 12.57 GeV, respectively 


Since a parton jet in the parton shower model is well 
defined, one can then calculate the double differential 
cross section for 7 — jet events. Shown in Fig. || are the 
normalized Ep distributions of the jet parton, 


1 dN' y—jet 

dEp 


( 11 ) 


with given E j, of the tagged direct photon, from HIJING 
simulations. As we can see that the transverse energy 
of the jet parton has a wide smearing around Ej due 
to the initial state radiations associated with the hard 
processes. Because of the rapidly decrease with Ep of the 
direct photon production cross section, the distribution 
is biased toward smaller Ep than Ej. The average Ep 
is thus smaller than Ej. In this paper, we will use the 
jet parton distribution in Eq. M to approximate the 
E^-smearing function in Eq. (§• 


If one triggers a direct photon with a given Ej,, one 
should average over the Ep smearing of the jet in the 
calculation of particle distributions in the opposite direc¬ 
tion of the tagged photon. Such a smearing is important 
especially for hadrons with pp comparable or larger than 

Ej. 

If we define the inclusive fragmentation function asso¬ 
ciated with a direct photon in pp collisions as, 


d Uz) = Y, t °( E t) [ dEp9„(Ep,Eq.) 


a,h 

K 1 

_T_ 7-)0 t * _ 

Ep h/a ^Ep' 


Ej, 


( 12 ) 


with z = pp/Ej, the hadrons’ momenta as fractions of 
the direct photon’s transverse energy, we can rewrite the 
pp spectrum [Eq. ([?])] of hadrons in the opposite direction 
of a tagged photon as 

dN;~ h± = DUpt/E^p) C{Ay,A<t>) 
dyd 2 pp PtE^ AyA<j> 

Considering parton energy loss in central AA collisions, 
we model the modified jet fragmentation functions as 
given by Eq. (j|). Including the Ep smearing and av¬ 
eraging over the 7 -jet production position in the trans¬ 
verse direction, the inclusive fragmentation function of a 
photon-tagged jet in central A + A collisions is, 



T, r ^ E r) J dE T g AA (Ep,El) 

' ' a,h 


Ep , El , 
x —D h/a (z—,L) , 

H/T H/t 


(14) 


where Taa ( 0 ) = J d 2 rt\(r) is the overlap function of AA 
collisions at zero impact-parameter. The Er-smearing 
function g AA (Ep, Ep) in A A collisions should be differ¬ 
ent from that in pp collisions due to initial multiple par- 
ton scatterings. However, for the moment, we will regard 
them as the same and postpone the discussion of the dif¬ 
ference to the next section. We have assumed that direct 
photon production rate is proportional to the number of 
binary nucleon-nucleon collisions. Neglecting expansion 
in the transverse direction, the total distance a parton 
produced at (r, <fi) will travel in the transverse direction 
is L(r,(j>) = y/R a — r 2 (1 — cos 2 </>) — rcos</>. Using the 
above inclusive fragmentation function in Eq. (13), one 
can similarly calculate the pp spectrum of particles in 
the opposite direction of a tagged photon in AA colli¬ 
sions. Setting the parton energy loss dE/dx = 0, the 
above equation should be reduced to the inclusive frag¬ 
mentation function in pp collisions in Eq. @ and the 
corresponding pp spectrum should also become the same 
as in pp collisions. 

To measure the modification of the inclusive fragmen¬ 
tation function in experiments, one should first select 
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events with a direct photon of energy E Then one 
measures the particle spectrum in the kinematical region 
(Ay, A <f)) in the opposite direction of the tagged photon. 
After subtracting the background which is essentially the 
Pt spectrum in ordinary events, one can use Eq. © to 
extract the inclusive jet fragmentation function, D 1 {z ), 
from the resultant spectrum. One can then compare the 
extracted inclusive jet fragmentation function in central 
AA collisions to that in pp or peripheral AA collisions to 
study the modification due to parton energy loss. Note 
that the centrality requirements for the signal (7 + jet) 
and background (ordinary) events should be the same. 
The overall acceptance factor C'(Ay, A<fi) in AA collisions 
remains approximately the same as in pp collisions with 
small but measurable corrections due to the kp broaden¬ 
ing of the leading parton as we will discuss later. 



p T (GeV/c) 

FIG. 3. The differential pp spectrum of charged particles 
from the fragmentation of a photon-tagged jet with =10, 
15 GeV and the underlying background in central Au + Au 
collisions at ifs = 200 GeV. The direct photon is restricted to 
\y\ < Ay/2 = 0.5. Charged particles are limited to the same 
rapidity range and in the opposite direction of the photon, 
|<() — — 7r| < A(f>/2 = 1.0. Solid lines are calculations using 

Eq. (i7|) and points are HIJING simulations of 20K events. 

To demonstrate the feasibility of the above prescribed 
procedure, we show in Fig. || the calculated pp distri¬ 
butions of charged hadrons from the fragmentation of 
photon-tagged jets with E^ = 10, 15 GeV and the un¬ 
derlying background from the rest of a central Au + Au 
collisions at the RHIC energy. We set dE/dx = 0 so 
the effect of parton energy loss is not included yet. The 
points are HIJING simulations of 20K events and solid 
lines for jet fragmentation are numerical results from 
Eq. © with the fragmentation functions given by the 
parametrization of e + e _ data ]l2| |. The numerical result 
(solid line) for the background coming from jet fragmen¬ 


tation in ordinary central events is obtained by the convo¬ 
lution of the jet cross section and fragmentation functions 
||. As we can see, the spectra from jet fragmentation 
are significantly higher than the background at moder¬ 
ately large transverse momenta. The background in pp 
collisions is about 1200 (the number of binary nucleon- 
nucleon collisions) times smaller than in central Au + Au 
collisions. One can therefore easily extract the inclusive 
fragmentation function from the experimental data with¬ 
out much statistical errors from the subtraction of the 
background. This conclusion remains valid even if one 
includes the parton energy loss in AA collisions because 
both the background and the particles from jet fragmen¬ 
tation are suppressed by approximately the same amount 
due to jet quenching [|]||. As one can expect from Fig. [|, 
for direct photons with E J < 6 GeV at the RHIC en¬ 
ergy, the hadron spectrum from the jet fragmentation is 
much smaller than the background. In this case, one can 
no longer accurately extract the effective fragmentation 
function with finite number of events. 

Shown in Fig are similar calculations for central 
Pb + Pb collisions at the LHC energy with Ej = 60 GeV 
for the tagged photons. It is clear that the overall back¬ 
ground is much larger than at the RHIC energy. There¬ 
fore, one needs to trigger on large E j, photons. Our calcu¬ 
lation shows that the fragmented hadron spectrum from 
the photon-tagged jets with Ej, = 40 GeV is roughly as 
large as the background at the LHC energy, from which 
one can barely extract the effective fragmentation func¬ 
tion. For Ej, < 40 GeV, the fragmentation spectrum is 
too small to be extracted. 

As a general criterion on the minimum value of E/ 
in our prescribed procedure in central AA collisions, one 
should require 

Taa ^^J^( E t™ m) = UGeV- 1 ). (15) 

Since large pp hadrons in the background also come 
from jet fragmentation, this is to ensure that Ej, is large 
enough such that the average number of jets with Ep = 
Ej, in each central collisions is less than 1. Only then, the 
inclusive fragmentation function of the photon-tagged jet 
can be extracted with confidence after the subtraction of 
the background. Shown in Table |, are values of E// min 
for different central A + A collisions at yfs = 200 GeV. 
This is consistent with what one can expect from Fig. [^. 


A 

80 

120 

160 

200 

&T min (GeV) 

5.0 

5.5 

6.0 

6.4 


TABLE I. The minimum transverse energy of the triggered 
photon, EJ min , required in order for the fragmentation func¬ 
tion of photon-tagged jets to be reliably extracted in central 
A + A collisions at y/s = 200 GeV. 
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p T (GeV/c) 

FIG. 4. The same as Fig. |d| except for Ef = 60 GeV at 
= 5.5 TeV. 

IV. EFFECTS OF JET E T SMEARING 

From Figs. || and [|, one also notices that there are sig¬ 
nificant number of particles with px larger than Ej, from 
fragmentation of the photon-tagged jets. This is because 
of the .Ex-smearing of the jet caused by initial state radia¬ 
tions. To illustrate the effect of the Et smearing, we plot 
in Fig. ^ the inclusive fragmentation functions (upper 
panel) with (solid lines) and without (dashed lines) Et 
smearing both for dE q /dx = 1 GeV/fm and dE q /dx = 0. 
The lower panel shows the ratios of the inclusive frag¬ 
mentation functions with and without energy loss. We 
assume that the mean-free-path of a quark is X q = 1 fm 
and the triggered photon has Ej = 15 GeV in central 
Au + Au collisions at y/s = 200 GeV. Notice that we now 
define 2 as a hadron’s fractional energy of the triggered 
photon. Because of the Ex-smearing of the jet caused 
by initial state radiations, hadrons can have px larger 
than Ej,. Therefore, the effective inclusive jet fragmen¬ 
tation function does not vanish at z = pt/E f >1. As 
we can see, the effect of Ex-smearing is only significant 
at large z. In particular at z ~ 1, the modified fragmen¬ 
tation function without Ex smearing has contributions 
only from those partons which escape the system with¬ 
out a single inelastic scattering, thus is controlled only 
by the mean-free-path Q. However, after taking into 
account of the Ex smearing, one also has contributions 
from the fragmentation of jets with Ex larger than E/f 
even if the jet has suffered energy loss. Therefore, the 
modification of the inclusive fragmentation function in 
this region of z depends on both the mean-free-path and 


the parton energy loss. Only at very large z > 2, the 
modification factor becomes independent of the energy 
loss, depending only on the mean-free-path. However, 
the production rate becomes also extremely small. For 
small and intermediate values of z, both the inclusive 
fragmentation function and the modification due to par- 
ton energy loss are not very sensitive to the Ex smearing. 



z 

FIG. 5. Upper panel:The inclusive fragmentation functions 
with (solid lines) and without (dashed lines) Et smearing for 
dE q /dx = 1 GeV/fm and dE q /dx = 0. Lower panel: The 
ratios of the inclusive fragmentation functions with and with¬ 
out energy loss. The mean-free-path of a quark is assumed to 
be X q = 1 fm and the triggered photon has Ef = 15 GeV in 
central Au + Au collisions at y/s = 200 GeV. 

To study the effect of jet Ex-smearing in detail, we 
show in Fig. ratios of the inclusive fragmentation 
function in central Au + Au collisions with energy loss 
dE q /dx = 1 GeV/fm over the one in pp collisions without 
energy loss. We shall refer to this ratio as the modifica¬ 
tion factor. The enhancement of soft particle production 
due to induced emissions is important only at very small 
fractional energy 2 . The fragmentation function is sup¬ 
pressed for large and intermediate 2 due to parton energy 
loss. For fixed dE/dx and 2 , the suppression becomes less 
as Ef or the average (Ex) increases. One can notice that 
there is an interesting structure in the region of z > 0.8 
which is also a consequence of the jet Ex-smearing. In 
this region, contributions from fragmentation of jets with 
Ex larger than Ef dominates. Since the leading parti¬ 
cles are relatively less suppressed for larger Et with fixed 











dE/dx and z, the decrease of the modification factor will 
then saturate at around z ~ 1 until the large Et tail of 
the Ey-smearing for the photon-tagged jet (as shown in 
Fig. ||) becomes insignificant. Therefore, the structure 
of the modification factor in the large z region results 
from the competition between the energy dependence of 
jet quenching and the falling-off of the tail of the jet Et 
smearing. The structure becomes more prominent for 
larger Ej, and it should also depend on the energy de¬ 
pendence of the parton energy loss. 

Another advantage of studying jet quenching in 
photon-tagged events is that the results are not sensitive 
to some of the effects of initial state multiple interac¬ 
tions, e.g., the nuclear shadowing of parton distributions 
(7J and energy loss of the beam partons, which can affect 
the direct photon production rate. However, there is one 
exception, i.e., the Et broadening from multiple initial 
and final state scatterings with the beam partons. Such 
Et broadening, which also causes the so-called Cronin 
effect [JUI in pA collisions, should also increase the Et 
smearing of the photon-tagged jets in both pA and AA 
collisions. This Et broadening has also been seen in dijet 
events of fixed target experiments |r| . Even though the 
Cronin effect for inclusive cross sections decreases with 
the colliding energy as indicated by current experiments 
[201 and theoretical estimates also predict it to be small 
[ L01 at the RHIC collider energy and beyond, the small 
Et broadening in the photon-tagged events should still 
have finite effects and one would like to have a handle on 
it experimentally. 



FIG. 6. The modification factor of the photon-tagged in¬ 
clusive jet fragmentation function in central Au+Au collisions 
at y/s = 200 GeV for a fixed dE q /dx = 1 GeV/fm 

One can directly measure the Et distributions of the 
photon-tagged jets in pp and pA collisions and study 
the possible change, using the conventional calorimetri- 
cal study of jets. However, due to small but finite back¬ 
ground in pp and pA collisions, it is difficult to mea¬ 
sure the calorimetric energy of jets with an accuracy 
of less than 1 GeV. Here we propose to study the Et 


broadening indirectly by measuring the modification fac¬ 
tor for the photon-tagged jets in pA collisions, since the 
inclusive fragmentation function is very sensitive to the 
Et smearing as we have demonstrated and the final jet 
partons in the central rapidity region do not experience 
transverse energy loss in pA collisions. In particular at 
2 = pt/E^, ~ 1, only those jets with Et > Ej con¬ 
tributes. The spectrum in this region is extremely sen¬ 
sitive to the Et broadening. Therefore, one should be 
able to measure even very small Et broadening via the 
modification factor in pA collisions. 

Let us assume that the transverse momentum kick 
from initial and final state scatterings with the beam par- 
tons has a Gaussian distribution with a width A pA , one 
can then calculate the effective jet fragmentation function 
similarly to Eq. © but with a modified Ey-smearing 
function g pA (Ey,Ej,), 


D; a (z) = ©r a (Ej) / dE T9pA(E T , 

a,h 


El) 


The modified E^-smearing function can be obtained as 
the convolution of the Ey-smearing function in pp colli¬ 
sions with a Gaussian distribution, 


gp A (E T , Ej.) = j dE a T J ^-e-^/ A ^ 5pp (E“,Ej) 

x S ^E t — \j Ef . 2 +Pt + 2ptEt cos 4 


dcj) 

tt 

2 / a 2 


e -Pr/ Ka 


e t r 7 n 2 

~7Ar~ 3 pp (Et j Ej) 

^pA 

Et 

\Je t 2 -Pt( 1 - cos 2 < j>) 


, (17) 


where E£ = y Et 2 — Pt (1 — cos 2 </>) — pt cos (j). For not 
very large values of A pA relative to Ej, the peak of the 
modified smearing function g pA (Er, Ej.) is simply shifted 
to larger values of Et as compared to < 7 pp (Et, EJ))) . One 
can characterize the E^-shift by 

AE T = J dE T E T [g p ^(E T , Ej) — g pp (E T , Ej)]. (18) 

To demonstrate the sensitivity of the effective frag¬ 
mentation function on the Et broadening due to multi¬ 
ple parton scatterings, we show in Fig [7] the modification 
factor Dj A (z) / D^ p (z) in pA collisions for three values of 
A Et with Ej = 10 and 15 GeV, respectively, at the 
RHIC energy. It is clear that the modification factor 
is sensitive to the additional Et smearing even for very 
small values of A Et- The shape of the modification fac¬ 
tor simply reflects the fact that the smearing function is 
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most modified around the peak Et ~ E J. Comparison 
between the calculation for E// = 10 and 15 GeV for fixed 
values of A Et shows that the relative effect of multiple 
scatterings decreases with increasing E J. 



z 

FIG. 7. The modification factor for the photon-tagged jet 
fragmentation function in pA collisions with E J = 10 and 15 
GeV at y/s = 200 GeV, for different values of A Et due to Et 
broadening. 

The Er-smearing function for AA collisions can be 
modeled the same way as in Eq. 0 ) for pA colli¬ 
sions, except that A p a is replaced by Aaa- According 
to the classical random-walk approximation |Q,^l|, A 2 a 
should be proportional to A 1 / 3 , or the average number 
of proton-nucleon subcollisions. In such an approxima¬ 
tion, A 2 aa = 2A? 4 . Using this modified Er-smearing 
function in Eq. (p"4|), we can calculate the modified effec¬ 
tive fragmentation function for the photon-tagged jets in 
AA collisions, including both the effect of parton energy 
loss through the dense medium and the additional Ep- 
smearing due to initial and final state scatterings with the 
beam partons. Shown in Fig. |] are the calculated modifi¬ 
cation factors with (solid line) and without (dashed line) 
Et broadening. It is clear that the Et broadening has 
significant effect on the final modification factor in AA 
collisions. One therefore has to study pp, pA and AA col¬ 
lisions systematically and subtract the effect caused by 
the Et broadening due to initial multiple scatterings to 
obtain the modification factor only due to parton energy 
loss in AA collisions. In the following discussions, we 
assume that such effect has already been subtracted out 


and we only concentrate on the effect of parton energy 
loss. 



0 0.5 1 1.5 

z=Pj/Ej 


FIG. 8. The modification factor for the inclusive fragmen¬ 
tation function of photon-tagged jets with (solid) and with¬ 
out (dashed) Et broadening due to initial parton scatter¬ 
ings, in central Au + Au collisions at y/s = 200 GeV. The 
parton energy loss is fixed at dE q /dx = 1 GeV/fm and the 
mean-free-path \ q = 1 fm. 

V. EXTRACTING PARTON ENERGY LOSS 

Given the modification of the inclusive fragmentation 
function of photon-tagged jets, one in principle should 
be able to extract the parton energy loss and the par- 
ton mean-free-path in our phenomenological model. The 
optimal case is when the average total energy loss is sig¬ 
nificant as compared to the initial jet energy, and yet 
the pt spectrum from jet fragmentation is still much 
larger than the underlying background. However, be¬ 
cause of the complication of the initial state radiations, 
one still cannot determine precisely the energy of the 
photon-tagged jet in each central AA event. Therefore, 
the parton energy loss and mean-free-path cannot be de¬ 
termined independently in a tangible way. As compared 
to our earlier results |J where we did not take into ac¬ 
count of the Et smearing of the photon-tagged jets, the 
modification of the averaged fragmentation function due 
to energy loss is quite sensitive to the value of the mean- 
free-path for dE q /dx = 1 GeV as shown in Fig. |(i|. 

To study the sensitivity of the modification to the en- 
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ergy loss, we plot in Fig. ^ the modification factor at a 
fixed value of z = 0.4 as functions of dE q /dx. For small 
values of dE q /dx , the suppression factor is more or less 
independent of the mean-free-path. This is referred to as 
the “soft emission” scenario in Ref. |j) where the suppres¬ 
sion is dominated by the leading parton with an average 
total energy loss (A Ej) = (n a )e a = ( L)dE a /dx . The 
suppression factor should scale with dE a /dx, depending 
very weakly on the mean-free-path. Assuming an expo¬ 
nential form of the fragmentation function ~ e~ cz for 
z = 0.2 ~ 0.8, one can show that the suppression factor 
has a form (1 — AE^,/E^,)exp(—czAE^/E^,). For large 
values of dE q /dx > 1 GeV/fm, the ratio is sensitive to 
the mean-free-path. However, as one can see from Fig. ||, 
the suppression factor flattens out as dE q /dx increases, 
especially for large values of the mean-free-path A a . This 
can be understood as the “hard emission” scenario in 
which the parton enegy loss per emission e a = X a dE a /dx 
is large. In this scenario, z n after n times emission be¬ 
comes so small, the n/0 contribution is completely sup¬ 
pressed in the modified fragmentation function in Eq. d). 
The only contribution is from n = 0 term, i.e., from the 
partons which escape the system without any induced 
radiation. In this case, the suppression factor is con¬ 
trolled by exp(—L/A q ), a factor independent of dE a /dx 
and Ej,. One thus needs to measure the suppression fac¬ 
tor at smaller values of z or a global fit to determine both 
the energy loss dE a /dx and the mean-free-path from the 
experimental data. 



dE q /dx (GeV/fm) 

FIG. 9. The modification factors for the inclusive 
fragmentation function of photon-tagged jets at given 
z = pr/E J = 0.4 as functions of parton energy loss dE q /dx 
with different values of the mean-free-path X q , in central 
Au + Au collisions at yfs = 200 GeV. 

Recent theoretical studies [Q of parton energy loss in a 
dense medium of a finite size L indicate that the energy 
loss per unit distance dE/dx could be proportional to 
the total distance that the parton has traveled since it is 
produced, 


dE a 

dx 


—T iln A 


-A, 


( 19 ) 


which was also shown to be proportional to the average 
transverse momentum broadening squared, A kj,, with re¬ 
spect to the direction of the initial parton momentum, 
where /i is the Debye mass of the medium and A the 
mean-free-path of the parton, C a =4/3 for quark and 
3 for a gluon. The kp broadening results from multiple 
scatterings which also induce the radiative energy loss for 
the propagating parton. 



FIG. 10. The modification factors for the inclusive 
fragmentation function of photon-tagged jets at given 
z = pr/E t = 0.4 as functions of A 1 / 3 in central A + A 
collisions at y/s = 200 GeV. The solid lines are for a dis¬ 
tance-dependent parton energy loss dE q /dx, while the dashed 
lines are for a constant dE q /dx. 

One way to test this experimentally is to study the 
modification factor at any given z value for different 
nucleus-nucleus collisions or for different centralities (im¬ 
pact parameters). Shown in Fig. are the modifica¬ 
tion factors for the inclusive fragmentation function at 
2 = 0.4 as functions of A 1 / 3 . We assume that the ra¬ 
dius of the cylindrical system is Ra = 1.2A 1 A fm. In 
one case (dashed lines), we assume a constant energy 
loss dE/dx= 0.5 GeV/fm. The modification factor de¬ 
creases almost linearly with A 1 / 3 . In another case (solid 
lines), we assume dE q /dx = 0.2(L/fm) GeV/fm. The 
average transverse distance a parton travels in a cylin- 
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drical system with transverse size Ra is (L) = 0.905-Ra- 
We choose the coefficient in dE q /dx such that its average 
value roughly equals to 0.5 GeV/fm for A = 20. To im¬ 
plement such an energy loss in our model, we assume the 
energy loss per scattering, for a parton traveling a total 
distance L , to be e a = A a 0.2(L/fm) GeV. As we can see, 
the suppression factor for a distance-dependent dE/dx 
decreases faster than the one with constant dE/dx. Un¬ 
fortunately, we have not found a unique way to extract 
the average total energy loss so that one could show 
that it is proportional to A 2 / 3 for the distance-dependent 
dE/dx. One possible procedure to determine the A de¬ 
pendence of the energy loss is to first determine dE/dx 
and A using a global fit to the measured modification fac¬ 
tor for each type of central AA collisions and then find 
the A dependence of the extracted dE/dx. However, such 
a procedure and our model depend on the assumption of 
the size of the dense medium produced in AA collisions. 

VI. Kt broadening and jet profile 

In our discussions so far, we have assumed that the 
jet profile in the opposite direction of the tagged photon 
remains the same in AA as in pp collisions, since we used 
the same acceptance factor C(Ay, Ac/). Such an accep¬ 
tance factor is determined by the effective jet profile in 
the opposite direction of the tagged photon. One can 
imagine that there should be two sources of corrections. 
One is due to the initial and final state multiple par- 
ton scatterings with the colliding nucleons. As we have 
discussed, such multiple scatterings can cause the broad¬ 
ening of the jet Et smearing. They shall also increase the 
acoplanarity of the jet with respect to the tagged pho¬ 
ton. One can study this effect directly via the effective 
jet profile in pA —► 7+jet + X processes as in dijet events 
& Let us assume that such increased acoplanarity can 
be measured and corrected. The second correction to the 
effective jet profile comes from multiple scatterings suf¬ 
fered by the leading parton while it propagates inside the 
dense medium. These multiple scatterings induce radia¬ 
tive energy loss and in the meantime also cause the kx 
broadening of the final parton with respect to its original 
transverse direction, giving rise to an additional acopla¬ 
narity . Such a change to the jet profile could affect the 
acceptance factor, which will be an overall factor to the 
measured jet fragmentation function if we assume the jet 
profile to be the same for particles with different frac¬ 
tional energies. 

Since we only consider jets in the central rapidity re¬ 
gion (y = 0), we assume that the final multiple scat¬ 
terings will only change the jet profile in the azimuthal 
direction. We define the jet profile function as f{(j)) = 
dEx/dip. If the initial effective jet profile is and 

the kx broadening distribution is given by a Gaussian 
form 0, the final effective jet profile function is then, 


POO 1 

f{4>) = J o dk 2 T ^-e- k T/ Ak Tf 0 (c/ - 4> jet ) , (20) 

where sin <j>j et = kx/Ex and A kj. is the average kx 
broadening squared. 



to the opposite direction of the tagged photon. The solid line 
is the original profile in pp collisions from HIJING simulations 
while the dashed line is the modified profile function with 
A k 2 T = 4 G eV 2 /c 2 . 

To demonstrate the effect of the kx broadening due to 
final multiple scatterings, we plot in Fig. [ll] (solid line) 
the azimuthal angle distribution of Ex (within \y\ < 0.5) 
with respect to the opposite direction of the tagged pho¬ 
ton with E j, = 10 GeV. We have subtracted the back¬ 
ground so that dEx/d(j) = 0 at (f> = n. The profile dis¬ 
tribution includes both the intrinsic distribution from jet 
fragmentation and the effect of initial state radiations. 
The acceptance factor is simply the fractional area within 
\<f>\ < A<j)/2 region. The kx broadening of jets due to mul¬ 
tiple scatterings will broaden the profile function. Shown 
as the dashed line is the profile function from Eq. ( |2(i| ) 
for Afcy = 4 (GeV/c) 2 . It is clear that with a modest 
value of the kx broadening, the acceptance factor only 
changes by a few percents. 

Since the change of the effective jet profile function is 
related to the average kx broadening, one can combine 
the measurement with the measured energy loss to verify 
the relationship between dE/dx and A kj. as in Eq. ©• 

VII. JET QUENCHING IN DEEP INELASTIC 
LEPTON-NUCLEUS SCATTERINGS 

Even though we have so far applied the parton energy 
loss in Eq. to a fast parton inside a dense matter, the 
generic form and its derivation is also valid for a parton 
propagating inside a cold nuclear matter or hot hadronic 
medium. The properties of the medium are manifested 
in the total transverse momentum kick Afc^. For a hot 
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QGP, Adirectly reflects the temperature, while for 
a cold nuclear matter it is related to the gluon density 
inside a nucleus ||. If there is a dramatic difference be¬ 
tween the transverse momentum broadening or the par- 
ton energy loss in QGP and a cold nuclear matter, then 
the measurement of parton energy loss in high-energy 
heavy-ion collisions can be used as a possible probe of 
QGP formation. It is thus also important to measure 
the parton energy loss in a cold nuclear matter. 

As we have mentioned in the Introduction, initial state 
interactions with beam nucleons prior to a hard pro¬ 
cess can also cause the participating partons to lose en¬ 
ergy, thus affecting the final cross section. Among many 
hard processes, such as Drell-Yan lepton pair and heavy 
quarkonium production at large x f in pA collisions [ p2| , 
the simplest processes where parton energy loss in cold 
nuclear matter can be directly measured are probably 
deeply inelastic lepton-nucleus scatterings. In such pro¬ 
cesses, one can relate the suppression of the leading 
hadrons to the attenuation of the quark jet inside the 
nuclear matter. There are many earlier studies of this 
problem in the literature |^3] ^b) . In this paper, we would 
like to revisit this problem within our framework of mod¬ 
ified fragmentation functions. 

In deeply inelastic IA collisions, a quark or anti-quark 
is knocked out of a nucleon by the virtual photon which 
carries energy v and virtuality Q. In the rest frame of the 
nucleus, the photon’s energy is transferred to the quark 
which then will propagate through the rest of the nucleus. 
If the hadronization time of the quark in the order of 
2f/Aq CD is much larger than the nuclear size, most of the 
leading hadrons from the jet fragmentation are formed 
outside of the nucleus. We can then attribute the atten¬ 
uation of the leading hadrons from the quark fragmen¬ 
tation to the energy loss of the propagating quark inside 
the nuclear matter. Let us assume that the longitudinal 
position of the nucleon from which the quark is knocked 
out is *|| in the direction of the virtual photon. Using 
Eq. (j|) for the modified fragmentation functions due to 
parton energy loss and averaging over the longitudinal 
and transverse position of the interaction point inside the 
nucleus, we can obtain the effective quark fragmentation 
functions in deeply inelastic lepton-nucleus collisions, 
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(z,Q 2 ) 
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^ Jo R\ 


r L A (x ± ) dx ^ 
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Dh/a(z,AL,Q 2 ) 


( 21 ) 


where La(x±) = \J R\ — x\ , A L = —xn + La(x _ l ), and 
a hard-sphere nuclear distribution is used. 

In a parton model, the nuclear structure function is 
defined as 


F% A (x, Q 2 ) = ZF* p/a {x, Q 2 ) + (A - Z)F* n/A (x, Q 2 ) , 
\F 2 £p/A (x,Q 2 ) = ^[2V n/a (x,Q 2 ) + 2 u n/a (x,Q 2 )) 


+ g [Vn/a(x, Q") + 2c?at/ j 4(*, Q~) + 2sn/a(x, Q~)\ , 
^F 2 n/A (x,Q 2 ) = ^[V n/a (x,Q 2 ) + 2u n/a (x,Q 2 )\ 

+ g[2VN/A( x ,Q 2 ) + 2d N /A(x,Q 2 ) + 2s N /A(x,Q 2 )] , ( 22 ) 

where Vn/a( x ,Q 2 ) (normalized to 1) is the effective va¬ 
lence quark distribution, qN/A( x ,Q 2 y s are the effective 
sea quark distributions per nucleon inside a nucleus and 
x = Q 2 /2rriNU. Here we neglect the isospin asymme¬ 
try in the sea quark distributions, i.e., un/a(x,Q 2 ) = 
d N /A(x,Q 2 ). Because of nuclear effects such as shadow¬ 
ing, the effective parton distributions per nucleon inside 
a nucleus are different from that inside a nucleon in the 
vacuum. There are many different mechanisms for the 
nuclear modification of the parton distributions and they 
could be different for valence and sea quark distributions 
|p7| . In this paper, we assume that the nuclear modifica¬ 
tion factors for quark distributions are the same and can 
be given by the ratio of structure functions, 


Ra/d{x,Q 2 ) 


Vn/a(x,Q 2 ) _ q N /A{x,Q 2 ) 
Vn/d(x,Q 2 ) cin/d{x,Q 2 ) 
2 F£ A (x,Q 2 ) 

AF 2 d (x, Q 2 ) 


(23) 



z 

FIG. 12. The suppression factor of charged hadron pro¬ 
duction in the jet fragmentation region in deeply inelastic IA 
collisions as a function of z = Eh/v. The lines are calcu¬ 
lations using modified fragmentation functions with parton 
energy loss dE q /dx. The data are from Ref. ps|. 


13 










With parton fragmentation model, one can also define 
the semi-inclusive structure function associated with pro¬ 
duction of hadrons with momentum zv in the direction 
of the virtual photon, 

F^ h± (x, z, Q 2 ) = ZF*»> ,A ^ h± (x, z, Q 2 ) 

+ (A-Z)F* n/A ^ h± (x,z,Q 2 ) , 


energy loss inside the nuclear matter, one still needs more 
accurate measurements in deeply inelastic £A collisions. 
Because of the finite total parton energy loss possible in¬ 
side a finite nucleus, events with small values of v (energy 
carried by the struck quark) are more desirable. 


^pip/A->h d 


— Fr 


\ln/'A—*hr 


(x, z, Q 2 ) = ^[2V N/A (x,Q 2 ) + 2u N/A (x,Q 2 )][D e n j v (z,Q 2 ) + D eA /v (z,Q 2 )\ 
+ \[V n/ a{x, Q 2 ) + 2 d N/A (x, Q 2 )][Dij v (z, Q 2 ) + D eA /s (z, Q 2 )} 
+ ^2 s N/A (x,Q 2 )[D eA s (z,Q 2 )+D eA /v (z,Q 2 )] , 

(x,z,Q 2 ) = ^[V n/a (x,Q 2 ) + 2u N/A (x,Q 2 )][D e n jy(z,Q 2 ) + D^ /v (z,Q 2 )\ 
+ i[ 2V n/a (x , Q 2 ) + 2d N/A (x, Q 2 )}[Di%{z, Q 2 ) + Dg /S (z, Q 2 )} 
+ ^2 s N/A (x,Q 2 )[Di%(z,Q 2 ) + D^ /v (z,Q 2 )\ , 


(24) 


In principle, one should also take into account hadron 
production from the fragmentation of the nuclear rem¬ 
nants. However, one can neglect them for relatively large 
values of z. In the above equation, we have used the fol¬ 
lowing definitions for the quark fragmentation functions, 


^tt/V — Dff/u D-x/u Dft/d D^/di 
Djr/S — D-n/s D-n/si 

Dk/v = D k / u = D k / u = D k / s = D k /s , 

Dk/s = D K/d = D K/d . (25) 


With the above model assumptions, one can study the 
quark energy loss by measuring the modification of the jet 
fragmentation functions via the ratio of the above semi- 
inclusive structure functions for different nucleus (A and 
-D-deuterium) targets, 


Ra/d{z) 


2F$ A ^ h± (x,z,Q 2 ) 

AF^ h± {x,z,Q 2 ) 


(26) 


Shown in Fig |l2], is the calculation of the above ratio 
within our model of the modified fragmentation func¬ 
tions, together with experimental data from E665 [p8| . 
In the experiments, the averaged values of Q 2 , v and x 
are different for different values of z. We have taken into 
account such kinematic effects in our calculation, espe¬ 
cially the nuclear modification of the quark distributions. 
In the experimental measurements, small values of z are 
correlated to small values of x, where nuclear shadow¬ 
ing of parton distributions is important. This is why the 
ratio Ra/d ( z ) is smaller than 1 even at small values of 
z. Within the errors, the data are consistent with our 
calculation with parton energy loss of dE/dx = 0.5-1 
GeV/frn. It is clear that in order to pin down the quark 


VIII. EXPERIMENTAL FEASIBILITIES 

To have an estimate of the experimental feasibility of 
our proposed 7 +jet measurement, we list in Table || the 
number of 7 + jet events per year per unit rapidity and 
unit (GeV) Et at the RHIC collider energy. We assume 
a central Au + Au cross section of 125 mb with impact- 
parameters b < 2 fm. We have taken a luminosity of C = 
2 x 10 26 cm _2 s _1 with 100 operation days per year. The 
7 + jet cross sections are taken from the compilation by 
the Hard Probes (HP) Collaboration j|ij As we can see, 
although the number of direct photons with E J = 7 GeV 
is large enough, the rate for Ej, = 15, 20 GeV is still too 
small to give any statistically significant measurement of 
the fragmentation function and its modification in AA 
collisions. If one can increase the luminosity by a factor 
of 10, the numbers of events for both E j, =10 and 15 GeV 
are significant enough for a reasonable determination of 
the fragmentation function of the photon-tagged jets. 

Given enough number of events, one still has to over¬ 
come the large background of 7r°’s to identify the direct 
photons. Plotted in Figs. [b| and 13 are the production 
rates of direct photons (solid line) and 7r°’s (dashed and 
dot-dashed lines) for central Au + Au collisions at the 
RHIC and LHC energies. The rate of ir° production is 
calculated with the same jet fragmentation functions em¬ 
ployed in this paper and convoluted with jet production 
cross sections ||. We can see that without jet quench¬ 
ing, 7r° production rate is about 20 times larger than the 
direct photons at pr = 10 GeV/c at yfs = 200 GeV. 
Fortunately, jet quenching due to parton energy loss can 
significantly reduce 7r° rate at large pt as shown by the 
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dot-dashed line. However, one still has to face 7 r°’s about 
3 times higher than the direct photons at px = 10 GeV/c. 
At larger px , the situation improves, but one loses the 
production rate. Since the isolation cut method normally 
employed in pp collisions to reduce the background for di¬ 
rect photons does not work any more, the only way one 
can identify them with high accuracies has to be through 
the means of improved detectors. 



FIG. 13. The spectrum of direct photon production (solid) 
as compared to 7r° spectrum with (dot-dashed) and without 
(dashed) parton energy loss ( dE q /dx = 1 GeV/fm, X q — 1 
fm) in central Au + Au collisions at y/s — 200 GeV. 


the LHC energy. We assume a luminosity of C = 2 x 10 2 ' 
cm^ 2 s ^ 1 with 50 operation days per year for Au + Au 
collisions. The production rates are reasonably high due 
to both the high luminosity and collider energy. However, 
the corresponding background of 7r°’s is also high (see 
Fig. @ which may make the detection of direct photons 
more difficult. 

At the LHC energy, y/s = 5.5 TeV, the production rate 


for Z° + jet becomes large even for reasonably large P/ ■ 
Listed in Table IV are the number of Z° + jet events per 
year per unit rapidity integrated over P'/ with differ¬ 
ent low cut-off values. The production cross sections are 
provided by T. Han based on calculations as described 
in Ref. Q. Note that the given Z° production rates 
are integrated ones, thus appearing to be larger than the 
differential rate of direct photon production at the same 
transverse momentum. One can detect Z° through the 
dilepton channel which has almost no background in the 
range of the dilepton invariant mass near Mzo. One can 
then apply the same procedure as we have discussed in 
this paper for direct photon events and measure the mod¬ 
ification of the effective jet fragmentation function due to 
parton energy loss. However, the drawback of using the 
dilepton channel of Z° decay is that the effective num¬ 
ber of events via this channel is about 6.7% of the total 
number of Z° events. 


(GeV) 

7 

10 

15 

20 

dN 1 ~ J et /dydEx /year 

20500 

3550 

400 

70 



TABLE II. Rate of direct photon production in central 
Au + Au collisions at y/s = 200 GeV, with luminosity 
C = 2 x 10 26 cnr 2 s _1 and 100 operation days per year. 


El (GeV) 

40 50 60 

dN 7 ~ J et /dydEr /year 

2880 1070 490 


TABLE III. Rate of direct photon production in cen¬ 
tral Au + Au collisions at y/s = 5.5 TeV, with luminosity 
£ = 2 x 10 27 cm~ 2 s _1 and 50 operation days per year. 


Pf (GeV) 

>20 >40 >60 

dN z - -' et /dy/year 

21100 7700 3470 


TABLE IV. Rate of Z° production in central Au + Au 
collisions at yfs = 5.5 TeV, with luminosity L = 2 x 10 2 ' 
cm" 2 s _1 and 50 operation days per year. 


Similarly, we also list in Table III the number of 7 +jet 
events per year per unit rapidity and unit (GeV) Et at 


15 


























IX. CONCLUSIONS 

In summary, we have studied systematically how one 
can measure the parton energy loss in 7 + jet events in 
central high-energy heavy-ion collisions within the frame¬ 
work of a modified fragmentation function model. We 
have demonstrated that an effective fragmentation func¬ 
tion of the photon-tagged jets can be extracted from the 
inclusive charged hadron spectrum in the opposite di¬ 
rection of the tagged photon. We have estimated the 
background from large pr hadron production to be small 
as compared to hadron production from the jet fragmen¬ 
tation for relatively large values of pr- We also provided 
estimates of the lower limits on Ej, in central A + A col¬ 
lisions in order for such extraction to be possible. We 
further show that the effective fragmentation function is 
sensitive to the parton energy loss possibly experienced 
by the parton during its propagation through the pro¬ 
duced dense matter. The sensitivity is characterized by 
the so-called modification factor via the comparison of 
the effective fragmentation function in AA with the one 
in pp collisions. 

We have explicitly taken into account the F/r-smearing 
of the photon-tagged jets for a fixed value of Ej, due to 
initial state radiations. We also demonstrated that the 
defined modification factors in pA collisions probe the Et 
broadening due to multiple initial and final state parton 
scatterings with the beam nucleons. One should subtract 
out the effect of such Et broadening when extracting 
the parton energy loss from the modification factor in 
AA collisions. We have also made detailed analysis of 
the modification factor within our model and studied the 
sensitivity to different forms of the parton energy loss, 
e.g., A dependence and the effective jet profile as a func¬ 
tion of the azimuthal angle in the transverse plane. We 
have also applied our model to jet quenching in deeply 
inelastic lepton-nucleus collisions from which one can ex¬ 
tract the parton energy loss inside a cold nuclear matter. 
Finally, we have examined the experimental feasibilities 
of our proposed study. Our estimates show that in order 
to have accurate measurements, one need somewhat in¬ 
creased luminosity about C = 2 x 10 2 ' cirU 2 s _1 at the 
RHIC energy. At the LHC energy, one can alternatively 
use Z° as the trigger and study the associated jet frag¬ 
mentation. 
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